Hydrodynamic lift of vesicles under shear flow in microgravity 
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The dynamics of a vesicle suspension in a shear flow between parallel plates has been investigated 
under microgravity conditions, where vesicles are only submitted to hydrodynamic effects such as lift 
forces due to the presence of walls and drag forces. The temporal evolution of the spatial distribution 
of the vesicles has been recorded thanks to digital holographic microscopy, during parabolic nights 
and under normal gravity conditions. The collected data demonstrates that vesicles are pushed 
away from the walls with a lift velocity proportional to jR 3 /z 2 where 7 is the shear rate, R the 
vesicle radius and z its distance from the wall. This scaling as well as the dependence of the lift 
velocity upon vesicle aspect ratio are consistent with theoretical predictions by Olla [J. Phys. II 
France 7, 1533-1540 (1997)]. 
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I. INTRODUCTION 

Vesicles are closed lipid membranes suspended in an 
aqueous solution. They can be considered as a simple 
mechanical model for living cells, especially simple ones 
like red blood cells. Their behaviour under shear flow has 
been the subject of several theoretical Jl|, [2|, y, |4| and ex- 
perimental [E El B El El studies. When the viscosity ratio 
between the inner and outer fluids is small [1], vesicles 
follow a tank treading dynamics in which the orientation 
of the main axis of the vesicle is constant and the mem- 
brane undergoes a tank treading motion. When they are 
sheared close to a wall, a lift force app ears which pushes 
vesicles away from the wall 0, [TO, Qj], [T2, EH • ^ gravity 
is perpendicular to the wall, vesicles reach a stationary 
altitude where the lift force is balanced by vesicle weight. 
Because of the rapid decrease of the lift force when the 
distance increases, the equilibrium distance under gravity 
is usually of order of the vesicle size @ . 

The unbinding and lift of vesicles from a wall have been 
investigated theoretically and numerically by several au- 
thors. Cantat and Misbah [ll| and Seifert [12[ studied 
the dynamical unbinding of adhering vesicles and the lift 
force when the distance between vesicles and the wall is 
small. Olla [10[ made a theoretical prediction of the drift 
velocity of a neutrally buoyant ellipsoidal vesicle when 
the distance to the wall is large compared to vesicle size. 
He found: 

d i = U(X,r 2 ,r 3 )^, (1) 

where 7 is the shear rate, R = (a^aa) 1 / 3 the vesicle's 
effective radius, ai, CL2, as the ellipsoid's semi-axes and 
z the distance between the vesicle center and the wall. 
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U is a dimensionless drift velocity which depends on the 
viscosity ratio A and the axis ratios r2 = ^2/^1 and rs = 
as/ai of the ellipsoid. Sukumaran and Seifert [T^ made 
a numerical study of the lift in the vicinity of a wall and 
found a similar form for the lift velocity. 

An experimental investigation of the unbinding and lift 
of vesicles under gravity was performed by Abkarian et 
al. [a]. By measuring the equilibrium height of vesicles 
of known density, they suggested the following empirical 
law for the lift force: 

R 3 

F z ~777— , (2) 

where h is the distance between the vesicle membrane 
and the wall (h c± z — R). Assuming a Stokes-like drag 
force Fd ~ r)(dz/dt)R, this would lead to a lift veloc- 
ity for neutrally buoyant vesicles that scales like jR 2 /h. 
This markedly differs from eq. [TJ However this experi- 
mental study is limited to short distances from the wall 
where lubrication effects are strong, and in a situation 
where gravity may also affect vesicle shape. This may 
result in very different scaling laws. To our knowledge 
no experimental study deals with the drift of vesicles at 
larger distances from a wall in a shear flow. 

The ESA project BIOMICS (BIOMImetic and Cellu- 
lar Systems) is dedicated to the investigation of the dy- 
namics of a suspension of vesicles under shear flow in 
microgravity conditions. Within the framework of this 
project several parabolic flight campaigns (ESA PF 43, 
CNES VP 59 and CNES VP 65) on board the Airbus 
A300-Zero G from Novespace (Merignac, France) were 
performed and an experiment will be part of the MASER 
11 sounding rocket payload. 

In this paper we present experiments, realized during 
the VP65 parabolic flight campaign, aiming at studying 
the lift of vesicles initially close to a wall under micro- 
gravity conditions where only hydrodynamic forces are 
present. A digital holographic microscope (DHM), work- 
ing in transmission with a partial spatial coherent source, 
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FIG. 1: Sketch and picture of the shear flow chamber devel- 
oped by SSC (color online). 

allows a 3D visualisation of the evolution of vesicle posi- 
tions in the field of view during microgravity 

The results show a drift of vesicles from the wall, at 
distances greater than vesicle size, which is compatible 
with lift forces given by eq. [TJ It underlines that the 
spatial distribution of vesicles across the channel thick- 
ness is strongly influenced by the shear flow and vesicle 
deformability. 

II. EXPERIMENTAL SETUP 

The experiments were performed in a Couette shear 
flow chamber designed and manufactured by the Swedish 
Space Corporation (SSC) for the MASER 11 sounding 
rocket flight foreseen for May 2008. As shown on Fig. [TJ 
it is made of two parallel glass discs, with a diameter of 
10 cm and a gap of 170 fim. The bottom disc is fixed 
and the top one can be rotated at constant speed by a 
computer controlled motor via a belt and a reducer, al- 
lowing shear rates between 0.5 and 50 s _1 . An inlet at 
the center of the bottom plate and several outlets at the 
periphery allow filling with experimental fluids, vesicle 
injection and rinsing. Vesicle suspensions and extra car- 
rier fluid are contained in syringes connected to the inlet 
via flexible tubes. A KDS syringe pump is used for the 
injections, and a soft pouch collects the excess and waste 
fluids. 

The observation system used to characterize the evo- 
lution and dynamics of the suspension is a digital holo- 
graphic microscope working with a partially coherent 
illumination designed and manufactured by MRC [141 ]. 
Digital holography is an emerging technology that pro- 
vides 3D imaging of microscopic samples. The main ad- 
vantage of digital holography is its capability to record 
rapidly the 3D information without scanning along the 
optical direction by changing the focus as with usual mi- 
croscopy or confocal microscopy. Therefore it provides 
the capability to investigate highly dynamic phenomena. 
This is a decisive advantage for 3D flow analysis. The 
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FIG. 2: Optical setup of the digital holographic microscope. 
LI: focusing lens; RGG: rotating ground glass for spatial co- 
herence reduction; L2: collimating lens; L3, L4: identical 
microscope lens (x 20); L5: refocusing lens; CCD: charge- 
coupled device camera; M1-M3: mirrors; BS1, BS2: beam 
splitters; S: shear flow chamber 



system takes benefit of a source of partial spatial coher- 
ence, created with a laser beam incident on a rotating 
ground glass that reduces the amount of coherent arti- 
fact noise and the multiple reflection perturbations [l5[. 
The optical system is described by the Fig. [2j It is 
a Mach-Zehnder interferometer in microscope configura- 
tion. The shear flow chamber is placed in one arm of 
the interferometer in front of a microscope lens. The 
direction of observation is perpendicular to the plates, 
approximately 3 cm from the axis of rotation. The mi- 
croscope lens, coupled with the refocusing lens, produces 
the image of one plane of the shear flow chamber thick- 
ness on the CCD. The second arm of the interferometer 
constitutes the reference beam and is also incident on the 
CCD where it is interfering with the object beam. 

By recording the holographic data resulting from the 
interference pattern, the digital holographic microscope 
provides the capability to refocus numerically vesicles 
that are out of focus with respect to the focus plane of the 
optical imaging subsystem. The holograms are recorded 
by a CCD camera and the digital holographic process, 
which consists of the refocusing in depth, is achieved 
by computation. The digital holographic reconstruction 
increases, without scanning, the depth of investigation 
by a factor of 100 in comparison with classical optical 
microscopy. The phase and intensity information are 
extracted from every hologram by the Fourier method 
[lrl Qj} and the complex amplitude is computed and used 
to refocus it in any plane inside the shear flow chamber 
thickness [18]. The field of view of the digital holographic 
microscope is (400 /im) 3 . The image acquisition system 
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Solution 


Content 


Viscosity 


Density 


interior 


10% glycerol 


1.51 cP 


1.045 g/1 




+ 200 mM sucrose 






exterior 


10% glycerol 


1.55 cP 


1.035 g/1 




+ 230 mM glucose 







TABLE I: Composition and physical properties of sample in- 
ternal and external solutions 



is able to record 24 holograms/second. As the vesicles 
are transparent objects, when they are focus, it is neces- 
sary to take benefit of the phase information, to locate 
them in the experimental volume. Digital holographic 
reconstruction coupled to specific algorithmic techniques 
for object detection [l8|, Qjl permits to obtain the size, 
projected shape in the XY plane, 3D position, and an 
estimation of the mean velocity of each flowing vesicle. 



III. SAMPLE PREPARATION AND 
EXPERIMENTAL PROCEDURES 

Vesicles were prepared by the electroformation 
method |20| , from lipids (1,2 - Dioleoyl-sn-glycero-3- 
phosphocholine (DOPC) from Sigma- Aldrich (P-6354)) 
deposited on ITO glass plates and an electroformation so- 
lution made of water, glycerol and sucrose (S-7903 from 
Sigma). Each electroformation provides a polydisperse 
sample with vesicle diameter ranging from 5 to 100 jam 
(median size between 15 and 25 /im), a volume fraction 
of a few percent and a total volume of about 2 ml. 

Samples were then diluted into a slightly hyperosmotic 
exterior solution made of water, glycerol and glucose (G- 
7528 from Sigma), which leads to a partial deflation of 
vesicles, a difference of refractive index necessary for visu- 
alization in the digital holographic microscope and a den- 
sity contrast responsible for sedimentation under gravity. 
The properties of internal and external solutions are sum- 
marized in table H The internal and external viscosities 
are nearly identical (viscosity ratio A = 0.99), and the 
density difference is 0.01 g/1, leading to sedimentation 
velocities under normal gravity of about 20 /xm/s for vesi- 
cles with a radius of 10 jam. 

The experimental procedure during flight is as follows: 
A parabolic flight consists in 31 parabolas, each providing 
22 seconds of microgravity (with a tolerance of ±5.10 -2 
g). Between parabolas (about 2.5 min), phases of normal 
and hyper gravity (1.5 to 1.8 g), during which the rota- 
tion of the shear flow chamber is stopped, allow a rapid 
sedimentation of the vesicles and provide a simple and 
reproducible initial condition with all vesicles lying on 
the bottom (fixed) disc. At the beginning of micrograv- 
ity phases, the rotation of the chamber starts with shear 
rates between 0.5 and 50 s _1 and images are recorded at 
24 holograms/s through the DHM. 




FIG. 3: Example of a phase map extracted from an hologram 
recorded in microgravity. White objects are vesicles, dark 
shapes are markers on glass plates (T-shapes on the bottom 
and V-shapes on the top plate). 



IV. DATA PROCESSING 

To obtain reliable statistics, algorithms and techniques 
have been developed in order to get automatic tools for 
the data processing. Indeed, for each parabola, 600 holo- 
grams are recorded leading to a huge amount of data. To 
obtain the 3D behaviour of vesicles flowing in the field of 
view, each hologram is reconstructed on the whole shear 
flow chamber thickness. From this set of reconstructed 
planes, the x and y position of each vesicle is obtained by 
phase contrast analysis. The z position (over the chamber 
thickness) is then determined using a specific metric [l8[ . 
Knowing the position of each detected object, their sizes 
and shapes (projected on the XY plane) are obtained by 
segmenting the phase map in the respective focus plane 
(Fig. i). 

Collecting the 3D behaviour of many vesicles passing 
in the field of view during microgravity time allows to 
gather statistical data about the evolution of the distri- 
bution of vesicles in the chamber thickness - due to lift 
forces and hydrodynamic interactions - as a function of 
time, as well as the influence of parameters such as vesi- 
cle size and deflation. Between 5 and 20 vesicles with a 
diameter greater than 10 /am are detected on each im- 
age. We get the spatial distribution and shape of several 
thousand vesicles for each microgravity period. 

Further data processing consists in computing each de- 
tected vesicle's effective radius R — (aia^ 1 / 3 ) /2 where 
a\ and <22 are the measured long and short axis of the 
projected vesicle shape (fitted by an ellipse) and its as- 
pect ratio ai/d2. Note that these are apparent (pro- 
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FIG. 4: Distribution of vesicle distances from the wall during 
microgravity for a shear rate of 50 s _1 , 5 < R < 50 /xm and 
1 < ai/a 2 < 2 

jected) parameters. A vesicle's real radius and aspect 
ratio are slightly bigger due to the inclination of tank- 
treading vesicles in shear flow p], 0, [To[ and can be re- 
covered by assuming a prolate ellipsoidal shape (where 
the two small axes are identical) and using the theoreti- 
cal relation between vesicle deflation and inclination from 
p], Qij. Then, filtering on size i?, aspect ratio axjai and 
averaging of results over a prescribed number of vesicles 
or consecutive frames can be performed. 

V. RESULTS AND DISCUSSION 

As the shear flow starts under microgravity, polydis- 
perse vesicles are lifted away from the bottom wall with 
drift velocities which depend on vesicle size and shape. 
An example is show on Fig. [4] where the vesicle cloud is 
seen to drift towards the shear flow chamber center. 

The large dispersion on this graph is due to several 
factors: the main factor is the polydispersity of the sam- 
ple (both in size and aspect ratio): in a typical sample, 
detected vesicle radii range from 3 to 50 jam and aspect 
ratios (long axis over short axis) between 1 and 2. This 
leads to very different lift velocities. Hydrodynamic in- 
teractions between vesicles also lead to a dispersion of ob- 
jects: as the suspension is sheared, vesicles move around 
each other and repel. Finally errors from the optical tech- 
nique and data processing can be estimated to ±3/im. 

For more quantitative results, it is therefore necessary 
to find a proper scaling and/or filter data to make results 
less sensitive to polydispersity, and to average results over 
several consecutive frames (or several detected objects 
with similar properties) in order to reduce the effect of 
noise in the data. 

Figure [5] shows the result of such a treatment for dif- 
ferent experiments at different shear rates, where only 
vesicles in narrow ranges of sizes and aspect ratios were 
kept. The effect of the shear rate is clearly visible: vesi- 
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FIG. 5: Average vesicle distance from the wall vs. time un- 
der microgravity conditions for different shear rates (filtering 
parameters: 5 < R < 10 /mi, 1.1 < a\/d2 < 1.25 ; averaging 
over 20 detections). 

cles drift much faster at higher shear rates 7. Note that 
no saturation is observed, indicating that vesicles are still 
quite far from the mid plane of the shear flow chamber 
(located at z = 85 /im) and that the lift force is still sig- 
nificant at a distance from the wall of order a few vesicle 
diameters. 

By integrating Olla's result (eq. [1]), one gets the fol- 
lowing scaling law: 

(|) 3 = 3t/7t (3) 

This suggests to redraw the results with rescaled pa- 
rameters (z/R) 3 and jt. An example is shown in Fig. [6] 
where a comparison of results in microgravity and normal 
gravity is made. Remarkably, in the absence of gravity, 
all results fall on the same straight line, in agreement 
with eq. [3l while under gravity a saturation appears, at 
a value which depends on shear rate. This corresponds 
to the case studied by Abkarian et al. [6[ where grav- 
ity balances the lift force as the latter decreases when 
the vesicle moves away from the wall. Note that signifi- 
cant scattering is still present on this figure, partly due 
to statistics and mainly to vesicle dispersion linked to 
hydrodynamic interactions. 

The slope of the data in Fig. [6] should be (3 = 3£7, a 
function of vesicle aspect ratio and viscosity. To deter- 
mine the dependence of f3 upon vesicle shape, Fig. [6] was 
redrawn by varying the range of aspect ratios used for 
the filtering of data, and a linear fit was performed. One 
expects vesicles with aspect ratios close to 1 (sphere) to 
have very small drift velocities (as the lift vanished for 
a perfect sphere in Stokes flow), while elongated vesi- 
cles should drift faster. The results are shown in Fig. [71 
where the slope /3 is plotted vs vesicle aspect ratio, and 
compared to Olla's prediction (with no fitting parame- 
ter) . The agreement is quite good both qualitatively and 
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bling transition. The study of the dispersion of vesicles 
in the z direction should bring interesting information 
about hydrodynamic interactions under shear. 

These different results show that after 20 s of micro- 
gravity, no equilibrium distribution seems to be reached, 
underlining the motivation for longer mierogravity times. 
As far as rheology is concerned, the observed depletion 
of vesicles near the wall should lead to a decrease of the 
effective viscosity of the whole sample, as this increasing 
vesicle-free layer acts as a lubrication layer where viscous 
dissipation is lower. 
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FIG. 7: Lift parameter f3 vs vesicle aspect ratio a±/a2 deter- 
mined from experimental data. Top horizontal axis: corre- 
sponding values of the reduced volume v. Dashed line: theo- 
retical prediction from Olla [l(| 
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FIG. 6: Rescaled vesicle- wall distance < (z 3 /R 3 ) > 
vs jt under mierogravity (top) and normal gravity (bot- 
tom) conditions for different shear rates (filtering parameters: 
5 < R < 50 /im, 1.1 < ai/<22 < 1.25 ; averaging over 15 
detections) . 

quantitatively in the range of available aspect ratios, ex- 
cept above 1.2 where the number of corresponding vesi- 
cles decreases and statistics are not as favorable. 

While Olla's model seems to work quite well for a vis- 
cosity ratio of 1 as shown in this paper, it will be interest- 
ing in the future to vary this parameter and see how the 
lift survives when approaching the tank-treading to tum- 
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